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Bulk Dynamics in Heavy Ion Collisions 
Peter A. Steinberg * 



O I The features of heavy ion coUisions that suggest the relevance of collective dynamics, 
as opposed to mere superpositions of nucleon-nucleon or even parton-parton collisions, 



are reviewed. The surprise of these studies is that bulk observables are far simpler than 
typical dynamical models of nucleus-nucleus collisions would imply. These features are 
shown to have a natural interpretation in terms of statistical-hydrodynamical models. The 
relevance of hydrodynamics to heavy ion collisions, coupled with the various similarities 
O i of the heavy ion data with that of more elementary collisions, raises very basic questions 



about its relevance to smaller systems. 



1. Introduction 



Recently, there has been dramatic progress in the understanding of the dynamics of 
heavy ion collisions. This is due to the availability of a large, high-quality data set 
^ ! spanning an enormous range in energy, rapidity and event centrality. With the turn- 
I on of the RHIC facility in 2000, we now have information about the basic features of 
these collisions up to yj/wv = 200 GeV, where ^/sN^ is the nucleon-nucleon center-of- 
mass energy. Although most of the analyzed data is from mid-rapidity (90 degrees in 
the center-of-mass system), all of the accelerator facilities have experiments dedicated to 
acquiring data over the full rapidity range (e.g. E895 at the AGS, NA49 at the SPS, 
and BRAHMS and PHOBOS at RHIC). Finally, the various heavy ion experiments are 
starting to converge on comparable centrality measures, making comparisons between the 
data sets as a function of the number of participating nucleons feasible. 

During the summer of 2004, the four RHIC experiments (BRAHMS, PHENIX, 
PHOBOS, and STAR) produced draft "white papers", summarizing their most important 
results and relevant interpretations [ d I2l El El ■ The result of these discussions indicates 
an ongoing paradigm shift both in the interpretation of the lattice results and the 
experimental data. The 20-25% deviation of the lattice energy density from the Stefan- 
Boltzmann Limit is no longer interpreted as approaching the weakly-interacting Quark 
Gluon Plasma. Rather, this deviation is presently understood (using N=4 SUSY QCD) 
as precisely the signature of a strongly-interacting plasma [ El E] • Thus, rather than 
expecting a weakly-coupled gas of quarks and gluons to be liberated in high energy 
heavy ion collisions, one might expect to create a strongly-interacting state instead, with 
completely different properties than previously expected. More importantly, one would 
not expect systems with these properties to be generated by the two body scatterings 
of asymptotic hadrons. The energy density achieved in these collisions imply particle 
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Figure 2. Total number of charged particles per 
Figure 1. dNch/ d7]\\^\^i/ {Npartf^) participant pair as a function of collision centrality 
compared with hadronic and fo^ three RHIC energies. d+Au data from = 
pQCD-based models. 2OO GeV is also shown. 



densities far too high for individual particle states to remain distinct. 

If this is the case, then one should expect the system created in heavy ion collisions 
to show collective behavior, characteristic of a system with zero mean free path, i.e. the 
hydrodynamic limit. This treatment of the system stands in stark contrast to the usual 
dynamical transport approach to heavy ion collisions. In this sort of approach, one expects 
multiple, independent stages, characterized by different dynamical mechanisms (shadowed 
parton distributions, parton production and reinteractions, quark recombination and 
chemical freezeout, fragmentation functions, hadron rescattering, thermal freezeout and 
hadronic decays). 

The independent combination of these various sources is encoded in models such as 
HIJING. There are also models that avoid the partonic stage altogether and proceed 
directly to the hadronic dynamics, such as RQMD. Interestingly, all of these models 
have difficulty reproducing the energy dependence of the particle density at mid-rapidity 
{dNch/dri\^r]\<i/ (Npartf^)) as shown in Fig. [T][[7j, which compares predictions for Au+Au 
and p + p from HIJING and RQMD, and compares them with data. While the models are 
tuned appropriately on p + p, both miss the Au+Au data systematically as a function of 
beam energy. Variables such as this integrate over all of the possible dynamical processes 
outlined above, thus showing the utility of bulk observables to be sensitive to both hard 
and soft processes, and thus early and late stages. Of course, a single observable at a 
particular energy and centrality is of little use. Rather, it is the combined systematics of 
energy, centrality, and rapidity which provide a handle on the most important dynamical 
contributions. 

The surprising conclusion, partially outlined in this proceedings, is that the systematics 
of bulk observables are far simpler than one might expect given the various mechanisms 
that could contribute to the final state. However, it seems that one must consider the 
features of particle production over the full phase space (and not just concentrate on 
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Figure 3. dN,h/dr]/ {Npart/2) in d+Au 
collisions at ^/snn = 200 GeV clS (X 
function of collision centrality. 



Figure 4. dN / drj/ {Npart/'^) for peripheral 
and central events at ^/snn = 19.6 and 200 
GeV. 



particular kinematic regions) to isolate the simple organizing features. 



2. Features of charged-particle multiplicities in A + A collisions 

The most striking feature of bulk particle production in heavy ion collisions is the 
approximate linearity of the total number of charged particles (Nch) with the number of 
participating nucleons {Npart)- Results onp + A collisions in the 1970's and 80's [|8j, found 
that Nch scaled linearly as Npart x N^. Recent PHOBOS data on d-l-Au collisions, shown 
in Fig. 121 [ El shows that this phenomenon persists at RHIC energies of ^/snn =200 
GeV. This simple behavior occurs despite non trivial changes to the particle density over 
the whole pseudorapidity range, shown in Fig. 01 What was not expected was that the 
same scaling would be present in Au+Au collisions from y^sJvTV = 20 to 200 GeV. This 
is a non-trivial result, in that the variation of dNch/drj with centrality, shown in Figs. El 
still maintains an overall constancy of {Net)/ {Npart/'^)- It should also be pointed out that 
while the d+Au data connects simply to the p + p data, the A + A data does not seem to 
extrapolate smoothly to the p + p limit. 

Another striking global feature is "limiting fragmentation" , the energy independence of 
the particle yields with energy when considering a system with a fixed collision geometry 
in a frame where one of the projectiles is at rest [ ^2] This feature was expected to be 
found in very limited regions of pseudorapidity, characteristic of the fragmentation of 
each projectile [E]- Instead, limiting fragmentation is observed over a relatively wide 
range in pseudorapidity, as shown in Fig. Moreover, one finds that the shape of the 
"limiting curve" depends only on the impact parameter, and once this is determined, 
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Figure 5. dN/drj' / {Npartf^) for peripheral 
and central events at three RHIC energies. 



Figure 6. Nch compared for A + A, 
e'^e~ and p + p collisions (both at ^/s 
and y/s/2. The data is compared to 
formulas described in the text. 



data from different energies lies on the same curve. Since the data peels off of the limiting 
curve at approximately the same distance from mid-rapidity, this phenomenon also seems 
to strongly constrain the mid-rapidity yields as well. Thus, the energy dependence of 
particle yields seems to be controlled by a global constraint that is obeyed at all available 
beam energies. 

The third striking feature of bulk particle production is the apparent universality of the 
total multiplicity for A + A, p + p and e~^e~ reactions at high energy, shown in Fig. IHl [ 
[121 Uni- Of course, this is only observed if one accounts for the "leading particle effect" 
in p + p collisions by working with an "effective energy" v/Sgg = y/s/2. The agreement of 
y4 + y4 and e^e~ without rescaling suggests on the contrary that there is no similar leading 
particle effect in heavy ion collisions. This is presumably connected to the fact that a 
typical participant is struck multiple times in A + A (unlike p + p and p + A collisions). 
The deviation of A + A and e+e~ at low energies can be understood as a consequence of 
the presence of a substantial baryochemical potential [|I3j, which tends to suppress the 
overall entropy (via the relationship S = {E — pV — ^bNb)/T. 

In conclusion, there seem to be three essential features of particle multiplicities in A-l-A 
collisions: 1) A'part-scaling of the total multiplicity, 2) A universal value of {Nch) I {Npart/2) 
in A + A, p + p and e~^e~ reactions, and 3) "limiting fragmentation" , which seems to 
constrain the global angular distributions and thus even the mid-rapidity density. Such 
simple behavior would not naturally be expected to arise out of a dynamical approach, 
where each stage is in-principle distinct from the others. For instance, it would seem 
fortuitous for nuclear shadowing, energy loss, and fragmentation functions to conspire to 
arrive at A^pa^rscaling from processes where semi-hard physics is controlled by the number 
of binary collisions, Ncou- It seems that any relevant physics scenario (such as the Color 
Glass Condensate described by L. McLerran [Ej) must capture these essential features. 
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Figure 7. BRAHMS data dN/dy for charged 
pions at 200 GeV, fit to a Gaussian. The 
inset shows the width divided by the Landau 
expectation as a function of ^/sNN. 



Figure 8. The Landau prediction 
for the multiphcity and rapidity dis- 
tributions, observed in the hmiting 
fragmentation reference frame. 



3. Hydrodynamic Descriptions of High Energy Collisions 

As discussed in the introduction, one of the major surprises from the RHIC data 
has been the relevance of relativistic hydrodynamics in the overall understanding of 
the experimental data. In some sense, this should not be so surprising. It is well- 
known that any system in which the mean free path approaches zero should show 
hydrodynamic behavior, even an apparently disparate system such as cold ^Li atoms [ I15j . 
This system can be initialized with a strongly asymmetric geometry and the subsequent 
pressure gradients lead to spatial deformations that have been successfully modeled with 
hydrodynamics. 

The relevant geometry in the early stages of a heavy ion collision (or even a. p + p 
collision) is one characterized by the nuclear radius in the transverse direction, and a 
large contraction in the longitudinal direction. Already one sees that the relevant time 
scales for longitudinal and transverse dynamics are very different, of order tt ~ R/cg 
(approximately several fm/c) in the transverse direction and tl ~ mR/{csy/s) (<^ 1 
fm/c) in the longitudinal direction. 

3.1. Landau Hydrodynamics 

Already in the early 1950's, Fermi and Landau considered a system with the geometry 
just described, replacing the potential complexities of a high energy nuclear colhsion 
with a slab of area 7tR\ and length A = RAvnj^fs (and thus volume V oc R^ / ^fsmq) 
in which all of the energy of the incoming projectiles is assumed to thermalize. This 
leads to an initial energy density ofeocs7VAr~4 TeV/fm^ at RHIC energies, a value 
typically seen as unphysically large. And yet, when one uses the massless blackbody 
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equation of state p = e/3 and converts this to a relation for the entropy density, a oc e^/^, 
the total entropy in the collision volume is = aV oc s^^^/^/s = s^/^, which is itself 
proportional to the total multiplicity It may not be immediately clear from 

these definitions, but in a nuclear coUision, S (y:V (x Npart, leading naturally to the Npart- 
scaling of total multiplicity. Landau extended this statistical model by using the equations 
of relativistic hydrodynamics [El CHI- It turned out that the strong compression along 
one axis leads to highly anisotropic distributions, with Gaussian rapidity distributions 
of width (J = In ^Jsjrr^ . It is the application of the Fermi-Landau initial conditions to 
the generally-accepted formalism of 3D relativistic hydrodynamics that is known as the 
"Landau hydrodynamical model", as advanced by Cooper et al. [EJ? Carruthers et al. [ 
OU] and Shuryak [|^, among others. What is striking is just how much existing data 
is broadly consistent with Landau's original predictions from the 1950's. As shown in 
Fig. El the Landau multiplicity formula, tuned on the low-energy data = 2.2s^/^ does 
a reasonable job on describing the trend of e~^e~ over a wide range of ^/s, and thus the 
p + p data at y/s/2 and the A + A data above y/s^N = 20 GeV [jUj. Of course it should 
not be overlooked that pQCD can describe the e"'"e~ data just as well, if not better, 
than the Landau formula. However, it is not clear why the two approaches agree to better 
than 10% over the range for which data exists, as shown in Fig. El 

The BRAHMS rapidity distributions of charged pions are distinctly Gaussian (shown 
in Fig. [7|) and have a width which deviates less than 10% from Landau's parameter-free 
predictions (as seen in the inset of Fig.[7j) It is less well-known that Landau's formulas 
{Nch oc s^/^ and cr = In ^JV/m?) actually predict limiting fragmentation when plotted as 
a function of y' = y — ybcam, as shown in Fig. |H1 [ISj- It appears that Landau's physical 
picture is already consistent with the essential features of particle multiplicities in heavy 
ion collisions. Moreover, in contrast to models that explain the rapidity distributions 
as the consequence of scatterings between the partons in the initial-state nucleon and 
nuclear wave functions [ 24, the Landau model starts with a static initial state and 
rapidly generates the rapidity distribution by means of hydrodynamics. Thus, dNch/drj is 
not a static "initial-state" effect, but rather the result of a dynamical process in the very 
first stages of the collision. 

3.2. Bjorken Hydrodynamics 

Although Landau hydrodynamics appears to be relevant for the physics in the very 
early stages (r <^ 1/m/c), it is still a model only understood semi-analytically with 
some particularly drastic approximations. By contrast, hydrodynamic calculations 
which assume boost-invariance in the initial conditions have been used for quantitative 
comparisons with experimental observables that are sensitive to early-time pressure 
gradients [(23- These models are based on the pioneering work of Bjorken [^7\, who 
postulated the imposition of boost-invariance [ |2E1 as a guiding principle for high energy 
reactions. Of course, since the calculations are initialized at time scales on the order of 
1 fm/c, they are unable to calculate the initial-state entropy. However, given this single 
piece of experimental data, and an assumed equation of state (usually a hybrid of the 
Landau EOS and a hadronic EOS, with a mixed phase), they are able to successfully 
calculate the effects of transverse pressure on particle spectra (radial flow) the mapping 
from the initial-state geometry to anisotropies in the final-state transverse momentum 
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Figure 9. f2 as a function of px for 
various particle species, compared to a 
hydrodynamics- inspired fit. 



Figure 10. f 2 as a function of rj compared 
with 3D hydro calculations. 



distributions (elliptic flow) [1^. 

Several clear signatures of radial flow are present in the experimental data even without 
recourse to particular models. At low transverse momentum [px < ^) it is observed 
that the particle spectra harden with increasing particle mass. This is characteristic 
of a collective flow velocity field that gives heavier particles a larger momentum kick 
(p = 7/3m). This immediately breaks the mr-scaling seen m p + p data, which is often 
thought to result from emission from a thermalized source. Models which include a more 
detailed handling of chemical freezeout are able to reproduce these trends [1^. 

Elliptic flow has now been comprehensively described by calculations using boost- 
invariant hydrodynamics, both as a function of centrality and as a function of pt and 
particle mass for a fixed centrality. An example is shown in Fig. |H1 which shows the 
characteristic "fine structure", or mass-splitting of the asymmetry parameter V2 as a 
function of the particle's transverse momentum [I2n]- These are non-trivial relationships 
between various species that are not typically described well in the dynamical approach 
described above. Nor can typical parton transport models reproduce the magnitude of 
the asymmetries [EDI- However, given the obvious discrepancy between the assumption of 
boost-invariance used in these calculations and the manifestly boost non-invariant particle 
distributions shown by the BRAHMS data, it is not surprising to find that it is difficult 
to reproduce the dependence of V2 on rj as measured by PHOBOS, shown in Fig. HU] [EI]- 
This is true even for truly 3D hydrodynamic calculations, showing the need to understand 
the initial conditions in some detail. The current state of the art calculations [ |221 use 
gluon structure functions calculated in the CGC framework, which explicitly avoids the 
longitudinal dynamics of the Landau approach. 

In conclusion, we have seen that the hydro approach appears to warranted by a wide 
range of data, although no existing model or code can describe every detail correctly. This 
is especially true when considering longitudinal dynamics, which has not yet been fully 
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Figure 11. dNchldrj' / {Npartf^) for A + A collisions compared to dN/dy'j, in e~^e reactions 
over a broad range of ^/s. 

incorporated into models that describe many features of the transverse dynamics. 

4. Relation to Elementary Systems 

Typically, the success of hydrodynamic models in heavy ion collisions is attributed to 
their being larger than hadronic length scales, which are thought to control the scattering 
processes by which the system equilibrates. And yet, it has already been shown that 
certain features of elementary collisions, such as the total entropy density, are manifestly 
similar to heavy ion collisions. Thus, it is not ruled out a priori to seriously consider the 
relevance of hydrodynamics to smaller systems, such as p + p and e^e~. 

The previous discussions have suggested that rapidity distributions are not merely 
the consequence of the nucleon or nuclear wave functions, but may be dynamically 
generated by Landau initial conditions. The overall similarity between the pseudorapidity 
distributions in A + A and e+e~ collisions at the same ^/s = 200 GeV (shown in Fig. ^2 
in the limiting fragmentation frame) certainly complicate separate explanations of these 
systems. It is also observed that the particle density per participant pair in A + A is similar 
to the density relative to the thrust axis in e~^e~ — * hadrons. Given these similarities, the 
observation of limiting fragmentation in p + p and e~^e~ (also shown in Fig. ^2 although 
less precise than that observed in A + A) is not surprising. 

However, the issue of whether or not these systems are truly equilibrated is contingent 
on whether or not collective behavior can be discerned. It is well-known that the particle 
yields in e~^e~ and p+p collisions are described by statistical models at the same chemical 
freezeout temperature as A + A at RHIC energies [ESI Ell- However, it is a subject of 
debate whether or not this is simply "phase space dominance" as opposed to a true 
equilibration process [EH]- Thus, the small systems are not regarded as truly collective. 
And yet, recent preliminary STAR data on identified spectra in p + p collisions, shown in 
Fig. 1121 [Eni, shows a qualitatively similar rise in the mean px as a function of particle 



9 




Figure 12. Preliminary data on {pj) 
measured m. A + A and p + p collisions 
at = 200 GeV. 



Figure 13. HBT radii as a function of kx 
for d+Au and Au+Au collisions divided 
by the radii measured in p + p collisions. 



mass. Other preliminary data on HBT correlations, shown in Fig. El [EZI, finds that 
the pronounced kx dependence of HBT radii in Au+Au, which is normally understood 
to yield information on radial expansion, is proportional to that found in d+Au and 
p + p. Thus, we are faced the situation that experimental observables are interpreted 
as expansion dynamics in the case oi A + A, but as merely reference data in the more 
elementary systems. 



5. Conclusions 

It is argued in this proceedings that soft physics in strongly interacting systems may well 
be simpler than typical dynamical models would generally suggest. In fact, hydrodynamics 
may provide a unified conceptual framework from the first moment of the collision all the 
way to freezeout. More than that, it may provide a basis for understanding elementary 
systems as well, perhaps complementing the approach of understanding these systems 
with perturbative QCD. Of course, QCD will be fundamental in understanding the basic 
degrees of freedom of these strongly-interacting systems and how they thermalize so 
quickly, as implied by the Landau approach. And yet, basic issues purely related to 
hydrodynamics, such as the integration of longitudinal and transverse dynamics will have 
to be tackled for real progress to be made. Finally, while the focus of RHIC physics has 
been at mid-rapidity, detailed understanding of the global features of the event across the 
whole rapidity range will be necessary to fully understand the collective nature of these 
collisions. 
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